Conductance quantization and transport gaps in disordered graphene nanoribbons 
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We study numerically the effects of edge and bulk disorder on the conductance of graphene nanoribbons. We 
compute the conductance suppression due to Anderson localization induced by edge scattering and find that 
even for weak edge roughness, conductance steps are suppressed and transport gaps are induced. These gaps are 
approximately inversely proportional to the nanoribbon width. On/off conductance ratios grow exponentially 
with the nanoribbon length. Our results impose severe limitations to the use of graphene in ballistic nanowires. 

PACS numbers: 73.23.-b, 73.50.-h, 81.05.Uw 



I. INTRODUCTION 



Measurements of electronic transport in graphene triggered 
an intense effort to understand the physical properties of 
this material.^ Both fundamental and applied aspects are 
currently being investigated by a large number of groups 
around the world. As the quality of the samples improves 
and other synthesis techniques are developed, the material 
changes from a regime where bulk disorder is the dominant 
electron-scattering mechanism at low temperatures to a bal- 
listic one, where boundary conditions, crystal alignment, and 
edge defects play a dominant role in setting the transport prop- 
erties. This regime is now experimentally accessible in ul- 
tra narrow ribbons, which are promising for developing high- 
frequency, low-noise, low-power field-effect transistors. 

Motivated by recent experiments 3 ^ and theoretical 
studies , 6 ' 7 ' 8 ' 9 ! 10 ' 1 1,12,13,14 in this paper we explore how edge 
roughness affects conductance in long nanoribbons with re- 
alistic widths, from several nanometers to tens of nanome- 
ters. Through numerical simulations we show that even very 
weak edge disorder has a marked effect in the conductance of 
these nanoribbons. Moderate amounts of edge roughness can 
substantially suppress the linear conductance near the charge 
neutrality point and induce a transport gap when the nanorib- 
bon is long and the number of propagating channels is small. 
This effect is a manifestation of quasi-one-dimensional An- 
derson localization. We compute transport gaps, localization 
lengths, and on/off ratios, and explore the combined effect of 
edge and bulk disorders. We also comment on the effect of 
inelastic scattering and dephasing. For nanoribbons with very 
weak edge roughness, conductance steps appear at noninteger 
values of the conductance quantum e 2 /h, irrespective of the 
lattice orientation. Our results indicate that producing quan- 
tum point contacts with current graphene nanoribbons will be 
extremely challenging and can only be achieved if either scat- 
tering at edges of the constriction is substantially suppressed 
or the edges themselves are defined with atomic precision. 
The paper is organized as follows. In Sec. |II]we describe the 
numerical method and present our results for the linear con- 
ductance of a graphene nanoribbon in the presence of edge 
disorder. In Sec. [HI] we show the results of our computations 



when bulk disorder is added to the nanoribbons and in Sec. [IV] 
we discuss the effects of weak edge roughness on the conduc- 
tance quantization steps of nearly ballistic nanoribbons. Con- 
clusions and final comments are left to Sec. [V] 



n. TRANSPORT IN THE PRESENCE OF EDGE 
DISORDER 

The simulations are based on the standard nearest-neighbor 
tight-binding model of a single-layer graphene, 2 The lin- 
ear conductance is evaluated through the recursive Green's- 
function technique. 15 An infinite nanoribbon is broken into 
three parts (see Fig. Q]): two (left and right) semi-infinite re- 
gions of width W modeling ideal contacts and a finite central 
region of length L where edge and bulk disorders are intro- 
duced. 
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FIG. 1 : (Color online) Schematic representation of the graphene 
nanoribbon setup used in the simulations. 

For a nanoribbon with perfect edges and no bulk disorder^ 
the conductance near the neutrality point Ep = can be zero 
(for semiconductor armchair) or a multiple of the conductance 
quantum 2e 2 /h (for metallic armchair and zigzag), depend- 
ing on the availability of states. The first discontinuity in the 
conductance appears when the Fermi energy Ep matches the 
minimum of first electron or hole subband; other steps are 
reached as the minima of consecutive subbands are crossed. 

In our simulations, edge defects are created by extracting 
lattice sites (carbon atoms) from both edges of the nanoribbon 
following a uniform probability distribution. It is assumed 
that atoms at the edges are always attached to two other car- 
bon atoms and passivated by a neutral chemical ligand, such 
as hydrogen. The control parameters are the number of etch- 
ing sweeps N swccp , which is related to the roughness ampli- 
tude, and the etching probability per site in the kth sweep, pk, 



which is related to the edge defect density. Unless otherwise 
specified, an averaging over ten realizations of each disorder 
case analyzed was carried out to decrease sample-to-sample 
fluctuations and facilitate the visualization of the results. 




FIG. 2: (Color online) (a) and (b): Energy dependence of the aver- 
age linear conductance of graphene nanoribbons with varying edge 
roughness. All nanoribbons have a the same length (L — 45 nm) 
and similar widths (W = 4.4 nm for armchair and W — 4.7 nm 
for zigzag), (c) Typical etching profiles used in (a) and (b) (only 
segments of the nanoribbon atomic structure are shown). Left: arm- 
chair and right: zigzag. 

In Figs. [2^ and [2J3 we show the linear conductance of 
nanoribbons with metallic armchair and zigzag lattice orienta- 
tions as a function of the Fermi energy for different values of 
7V sweep and {pk} at zero temperature. Also shown in Fig. [2J: 
are the typical etching profiles for the four different cases. We 
have defined the roughness parameter r = (W — W)P/ao, 
where W is the average ribbon width, P = ~^2 k Pk, and is 
the lattice constant. 17 While a staircase of conductance steps 
is seen in the absence of edge roughness, the conductance 
rapidly degrades as the concentration and depth of the random 
etching increase. The conductance is strongly suppressed near 
the neutrality point even for relatively shallow etchings. Close 
to the neutrality point, a deep gap develops. 

The formation of a transport gap in ultranarrow graphene 
nanoribbons is very much consistent with all experimental ev- 
idence available so fa n 3,18 ' 19 Several mechanisms have been 
proposed to account for this phenomenon, such as straight- 
forward lateral confinement to many-body effects. 20 ' 21 Some 
of these mechanisms require the existence of substantial edge 
disorder along the ribbon (enough, for instance, to form bottle 
necks and quantum dots. 21 ) 




FIG. 3: (Color online) Average conductance of (a) metallic arm- 
chair (W = 4.3 nm) and (b) zigzag (W = 4.6 nm) nanoribbons 
as a function of the Fermi energy for three different widths when 
moderate roughness is present (zero temperature). The dashed lines 
indicate the conductance gap estimated by the change in the slope 
of the curves (see insets). A total of ten realizations of the edge 
roughness types #3 for armchair and #4 for zigzag were used for 
each curve presented. 



The results presented in Figs. an d Et show that the 
suppression of conductance can also occur at small values 
of r. Near the neutrality point, the number of propagating 
channels in the nanoribbon is very small and the system be- 
haves as a quasi-one dimensional wire. Edge defects act as 
randomly positioned short-range scatterers and induce strong 
backscattering, which in turn leads to Anderson localization 
if the nanoribbon is longer than the localization length. Note 
that even at room temperatures, we expect dephasing lengths 
in graphene to be exceedingly long. 22 Therefore, in practice, 
localization lengths can be shorter than both the nanoribbon 
length and the dephasing length. 

Further evidence of this effect is provided in Figs. [3^ and 
I3J) where the conductance as a function of Fermi energy is 
shown for three different nanoribbon widths. At a fixed en- 
ergy, for nanoribbons with increasing width, the localization 
becomes weaker and, consequently, the conductance suppres- 
sion decreases. In order to make contact with experiments, we 
estimated transport gaps by determining the energy value at 
which the curves show an inflection point [see insets of Figs. 
[3^ and03]. This inflection point corresponds to the crossing 
between two straight lines: an energy-independent conduc- 
tance at low energies and a linear function at high energies. 
The results are shown in Fig. @h. For both lattice orientations, 
we find that the transport gap E g scales approximately with 
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FIG. 4: (a): Transport gap dependence on the nanoribbon width. 
The dashed and dashed-dotted lines are fittings to the data as ex- 
plained in the main text, (b) On/off linear conductance ratio as a 
function of ribbon length. The solid lines are guides for the eyes 
and the dashed lines indicate fittings to exponential curves. Each 
data point corresponds to an average over 50 realizations of edge 
roughness. Inset: localization length as a function of ribbon width 
(circles are for armchair and squares are for zigzag). Solid (empty) 
symbols represent Ef ~ (Ef = 0.2 1). Lines are fittings of the 
functional form £ = A W a (see main text). 



FIG. 5: (Color online) The effect of combined edge and bulk dis- 
orders on the average conductance of short graphene nanoribbons. 
(a) Armchair and (b) zigzag orientations (both with the same width 
W = 8.3 nm and length L = 22 nm). The edge disorder corre- 
sponds to the case #2 of Fig. [2};; the bulk disorder has parameters 
fiimp = 0.04, Kq — 0.5, and £ = 2 ao. Inset: average linear con- 
ductance of armchair nanoribbons with bulk disorder and perfect 
edges at zero temperature (L = 22 nm, W = 4.4 nm, and five 
realizations per curve). The dashed lines correspond to correlation 
lengths £/ao = 10, 3, and 1 from top to bottom (ni mp = 0.02 and 
Kq — 1 in all three cases). 



the inverse of the nanoribbon width W and is only weakly de- 
pendent on the length L provided the latter is sufficiently long 
(not shown). Notice that the value obtained for A, the scaling 
prefactor, is in the same range of those found experimentally 
(A w 0.2 - 0.6 eV-nm) (Refs. OJandOl) even for moderate 
roughness. The gap is less pronounced for zigzag nanorib- 
bons, since for this orientation most of the current at low dop- 
ing is carried through bulk states 6 which are less sensitive to 
edge defects. Nevertheless, once the etching goes deeper than 
one lattice spacing, a clear gap develops for zigzag orienta- 
tions as well. 

In Fig. |4j} we plotted the ratio between the on and off lin- 
ear conductances as a function of the nanoribbon length. The 
off conductance G s was obtained at Ep = while G on 
was defined as the conductance at the Fermi energy where a 
transition from the first to the second step occurs in the clean 
nanoribbon. The curves indicate an approximate exponential 
growth in the on/off ratio, consistent with the Anderson lo- 
calization picture. Since conductance is broadly distributed 
in the localized regime, the on/off ratio develops very large 
fluctuations when the nanoribbon is long. In the inset we 
show the localization length l\ oc extracted by fitting to the 
data an expression of the form G(L) — G(0) e~ L / iloc . We 
note that the localization length grows with increasing ribbon 
width and energy. Fittings to the form Zi oc = AW a yield 
a = 0.21 - 0.61 at E F = and a = 1.8 - 2.2 at E F = 0.2 1, 
with the prefactor in the range A = 0.04 — 7.3 for Zi oc and W 
in lattice units. Thus, l\ oc can be comparable to W. 



ing along the nanoribbon N- imp Gaussian scatterers of width £ 
with random amplitudes { U n } drawn from a uniform distribu- 
tion [-SU, <5/7]. The intensity of the disorder is characterized 
by the dimensionless parameter Kq, which is defined through 
the correlation function 



(J/imp(rt)E4n P (r 3 -)) = K 



hv 



-|r<- 



2 /2? 2 



(1) 



In the dilute limit, when Ni mp is much smaller than the to- 
tal number of sites in the nanoribbon, one can show that 24 
Kq w 40.5ni mp (<5/y/£) 2 (c;/ao) 4 , where ni mp is the scatterer 
density per lattice site. For large graphene sheets at high 
doping, far from the neutrality point, it is possible to relate 
this parameter to the transport mean-free path using the Born 
approximation: 25 £ tI = 2Xf/(ttKq), where Xp is the Fermi 
wavelength in the graphene sheet, with <C £tr- 

The difference between armchair and zigzag states near the 
neutrality point also explains the effect of bulk disorder in 
nanoribbons with rough edges (Fig. |5j. For the metallic arm- 
chair orientation the low-lying states are concentrated at the 
edges and are quite sensitive to edge roughness. 6 In this case, 
moderate bulk disorder has a small effect on the transport gap. 
For zigzag orientations, the situation is the opposite, as bulk 
disorder disrupts the current-carrying states and suppresses 
conductance in this case. However, for both orientations, bulk 
disorder only leads to strong localization when £ < ao (short- 
range disorder). 



III. BULK DISORDER 



IV. CONDUCTANCE QUANTIZATION 



In order to investigate the effects of bulk disorder, we added 
an on-site correlated Gaussian disordered potential [/; mp (r) to 
the nanoribbon. 23 ' 24 ' 25 The latter is constructed by distribut- 



While it is clear that moderate edge disorder leads to sub- 
stantial suppression of the conductance of nanoribbons, what 
happens when the etching is nearly perfect and only very few 



4 




FIG. 6: (Color online) Suppression of conductance quantization 
steps due to weak edge disorder, (a) Armchair and (b) zigzag ori- 
entations. A total of 100 realizations were used in each case. Solid 
(dashed) lines indicate zero (room) temperature in the contacts. 
Inset: average conductance versus energy for a metallic armchair 
nanoribbon with L — 22 nm, no bulk disorder, p\ — 0.05 (ten 
realizations), and increasing widths: W = 4.4, 9.1, and 18.5 nm; 
solid (dashed) lines represent zero (room) temperature, while dotted 
lines describe the perfect edge case at zero temperature. 



and shallow defects are present? How does weak edge rough- 
ness affect conductance quantization in comparison to bulk 
disorder? Answers to these questions are provided in Fig. [6] 
and in the inset of Fig. [5^. The main effect of very weak edge 
disorder is to lower the conductance steps without changing 
their positions in energy. Some fluctuations occur near the 
transition regions between steps because of the sensitivity of 
the evanescent modes to variations in nanoribbon width, but 
these fluctuations are washed out at finite temperatures. How- 
ever, thermal fluctuations, even at room temperature, do not 
change the wider steps, as can be seen in the insets of Fig. [6j 
as long as the nanoribbon is sufficiently narrow. 

Bulk disorder, on the other hand, has a quite distinct ef- 
fect on the conductance quantization. In the inset of Fig. 
13 a) we show how bulk disorder affects the conductance steps 
of a metallic armchair nanoribbon. As the disorder range 
widens, the steps are smeared without shifting the conduc- 
tance value. This result can be understood in the following 
way. When only weak edge disorder is present, the width of 
the nanoribbon is hardly unaltered and propagating channels 
in the nanoribbon open up at the same energies as in the case 
of perfect edges. Yet, backscattering due to randomly posi- 
tioned edge defects, albeit weak, reduces the overall conduc- 
tance and shifts the steps downward. Long-range bulk dis- 
order, on the other hand, creates potential inhomogeneities 
which lead to the appearance of electron and hole puddles 
when the Fermi energy lies close to the Dirac point. 23 Trans- 
mission through these puddles creates mode mixing, which in 
turn smears the conductance steps. 

A simple model can be used to describe the suppression of 
the conductance steps. For the first step, let us assume that 
carriers propagate in one dimension through a sequence of 
randomly positioned but identical barriers. When the barrier 
reflectance is very small, R <C 1, it is straightforward to show 
that, in the short-wave limit (A <C L), the conductance in the 
first step goes as G\ ~ (2e 2 /h) [1 — NR], where N is the 



number of barriers, which can be directly related to the de- 
fect probability or density as follows: N = (L/a^)pi. For 
higher steps, more than one propagating mode is present and 
the system becomes quasi-one dimensional. In the absence 
of mode mixing, the conductance of the nth step behaves as 
G n w n (2e 2 /h) [1 — NR]. However, mode mixing is un- 
avoidable for n > 1 and one expects strong deviations from 
this simple scaling behavior 

The suppression of the conductance step as a function of 
nanoribbon length and edge defect concentration is presented 
in Fig. [7] For the first step of the armchair orientation the sim- 
ple one-dimensional scattering model works quite well. Both 
dashed lines in Figs. [7^ and |7p correspond to R = 0.035. 
The linear scaling ceases to apply for higher steps, with a sub- 
linear dependence indicating substantial mode mixing. For 
zigzag orientations, the suppression of the first conductance 
step was too small to be shown on the same plot. Figures [7J; 
and|7}i show that the simple scaling behavior no longer applies 
already when n = 2 (R = 0.013 was used in this case). 




FIG. 7: (Color online) Relative suppression of the nth conductance 
step as a function of (a) and (c) nanoribbon length and (b) and (d) 
defect concentration for metallic armchair (W = 4.4 nm) and 
zigzag (W = 5.2 nm) orientations. Each data set corresponds to 
an average over 100 realizations. The index n indicates the order of 
the conductance step and the solid lines are guides for the eyes. The 
dotted lines are explained in the main text. 



V. CONCLUSIONS 

Our results indicate that creating graphene quantum 
point contacts will depend fundamentally on atomic scale 
engineering. 26 In searching for ways to improve the nanorib- 
bon conductance in the presence of edge roughness, we found 
that, due to the Klein tunneling effect, side gates are not effec- 
tive in reducing edge scattering (i.e., electrostatic potentials do 
not confine Dirac fermions). At moderate roughness, Ander- 
son localization develops and a transport gaps appear. There 
is a quantitative agreement between the gaps that we find nu- 
merically when extrapolating our results to wider nanoribbons 
and the available experimental data, although we expect that 
other mechanisms, such as spin gaps 20 and charging effects 21 
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(not taken into account in our calculations), will likely com- 
pete with localization. Electronic dephasing can be introduced 
in the calculation following a standard procedure. 2 ? Since, in 
practice, dephasing lengths exceed the nanoribbon width, we 
expect that the main effect of dephasing will be the appearance 
of an additional, weakly energy dependent, suppression of the 
conductance, as the nanoribbon will break into a series of in- 
dependent quantum resistors. At this point, we would also 
like to note that we recently became aware of similar work by 
Evaldsson et alS^ 
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